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Interference represents one of the most striking manifestation of quantum physics in low¬ 
dimensional systems. Despite evidences of quantum interference in charge transport have been 
known for a long time, only recently signatures of interference induced thermal properties have been 
reported, paving the way for the phase-coherent manipulation of heat in mesoscopic devices. In 
this work we show that anomalous thermoelectric properties and efficient heat rectification can be 
achieved by exploiting the phase-coherent edge states of quantum Hall systems. By considering a 
tunneling geometry with multiple quantum point contacts, we demonstrate that the interference 
paths effectively break the electron-hole symmetry, allowing for a thermoelectric charge current 
flowing either from hot to cold or viceversa, depending on the details of the tunnel junction. Cor¬ 
respondingly, an interference induced heat current is predicted, and we are able to explain these 
results in terms of an intuitive physical picture. Moreover, we show that heat rectification can be 
achieved by coupling two quantum Hall systems with different hlling factors, and that this effect 
can be enhanced by exploiting the interference properties of the tunnel junction. 


I. INTRODUCTION 

In recent years a great attention has been devoted 
to the study of thermal transport at the nanoscale 
and energy and heat exchanges in small quantum de¬ 
vices m\- A deep understanding of these phenomena is 
of paramount importance for applications in solid-state 
cooling, high-precision sensors, cryogenic thermometry 
and thermal logic in quantum information 011 uni E]. 
This emerging field goes under the name of caloritron- 
ics Emm. Here, surprising experimental results have 
been already obtained. Among them Giazotto et al. have 
recently demonstrated the possibility to coherently ma¬ 
nipulate the heat flux in a hybrid superconducting circuit 
By realizing the thermal version of the electric 
Josephson interferometer, they have paved the way to¬ 
wards the phase-coherent manipulation of heat in meso¬ 
scopic devices m- Moreover, implementations of meso¬ 
scopic heat engines and thermal diodes have been pro¬ 
posed and will be soon achieved Ennui. These can 
be realized using quantum dots [H [2Qll22] , optomechan¬ 
ical systems m, and multiterminal mesoscopic devices 
combining also normal metals, superconductors and fer- 
romagnets [24U26] . The rapid progresses made in the 
field of caloritronics allow to glimpse a future in which 
electric and thermal manipulation will proceed on equal 
footing. 

In this context a promising role is played by topologically 
protected states, such as the edge states of quantum Hall 
systems (QHS) and topological insulators (TI). As far as 
the recently discovered TIs are concerned, the presence of 
protected helical edge states not only allows to generate 
peculiar spin-dependent thermal phenomena [271 EH] , but 
also to achieve high thermoelectric performances [29ll32] . 
On the other hand, the interest of the scientific commu¬ 
nity for QHS has been refueied in the view of possibie 
thermai appiications [71[24l[33l[34|. For exampie, Sanchez 


et al. have demonstrated that a three terminai device in 
the quantum Haii regime can work as a perfect thermai 
diode, with a rectification coefficient rg ^ oc, expioiting 
the chiraiity of edge states m- Heat transport measure¬ 
ments were also proposed in order to extract important 
information on fractional statistics and neutral modes in 
exotic fractional quantum Hall states |33l |34|. In these 
topologically states protection from backscattering guar¬ 
antees phase-coherent ballistic transport of charge and 
heat over long (> jam) distances [35H39] . and the possi¬ 
bility to control their charge transport properties in the 
presence of applied bias voltages via interferometric se¬ 
tups have been extensively studied. However, the effects 
of quantum interference on their thermal and thermoelec¬ 
tric properties remain quite unexplored, and our work, 
motivated indeed by the recent interest in phase-coherent 
heat manipulation, moves in this direction. 

In this paper we show that the chiral edge states of 
QHS can be exploited to implement coherent caloritron¬ 
ics. In particular, we consider two separated QHS cou¬ 
pled by a tunneling region, driven out of equilibrium by 
the presence of a thermal gradient that induces finite 
charge and heat flows. We demonstrate that a tunnel 
junction with n quantum point contacts (QPGs) enables 
to control the charge and heat transfer between the two 
QHS. As far as the charge sector is concerned, we show 
that, by varying the geometrical parameters of the junc¬ 
tion, particle-hole symmetry can be broken, allowing to 
selectively enhance electron tunneling with respect to 
hole tunneling or viceversa. Therefore we predict the 
interesting result that quantum interference can be ex¬ 
ploited to selectively switch the charge flow induced by 
the thermal gradient, that is, charge current can either 
flows from hot to cold or from cold to hot, depending on 
the interference properties. 

Remarkably, we find that quantum interference phenom¬ 
ena control heat transport as well. Contrary to the 
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charge transport, with electrons and holes giving oppo¬ 
site contributions, heat transport does not depend on 
the charge of the carriers, so that electrons and holes 
equally contribute; this property manifests in the corre¬ 
spondence between the zeros of the charge current and 
the maximum/minimum visibility of the heat current. 
Furthermore, we study how these quantum interference 
phenomena, already present in the absence of interac¬ 
tions (tunneling between integer QHS), are affected by 
the presence of interactions, by considering tunneling of 
electrons between fractional and integer QHS, where e-e 
interactions play an important role. As a general remark, 
a strong suppression of the signal appears, due to the 
anomalous temperature dependence of the effective tun¬ 
neling density of states, a hallmark of non-Fermi liquid 
behavior. More interestingly, the presence of different 
filling factors breaks the left-right symmetry, inducing 
rectification effects. By taking advantage from the inter¬ 
ference patterns induced in the presence of several QPCs, 
we show that strong rectification effects can be obtained. 
We thus demonstrate that the interplay between inter¬ 
actions and quantum interference is crucial in order to 
enhance the heat rectification. 

The paper is organized as follows. In Sec. |TI| we de¬ 
scribe the setup and evaluate the charge and heat cur¬ 
rents for a generic tunneling region. Sec. m is devoted 
to th e main discussion, focusing on in terfere nce phenom¬ 
ena (IIIA) and rectification effects (IIIB) in a multi¬ 
ple QPC geometry. Finally, we draw our conclusions in 
Sec.IVl 


II. MODEL AND TRANSPORT EQUATIONS 



Figure 1. (Color online) Scheme of two fractional quantum 
Hall systems with different filling factors ul and ur at temper¬ 
atures Tl > Tr. Counterpropagating edge states are coupled 
by a tunneling region, schematically depicted in the middle. 
The zoom in the right panel shows the case of tunneling due to 
multiple quantum point contacts (n = 3 in this case) equally 
spaced with distance d. 

We consider two quantum Hall bars with filling fac¬ 
tor {a = R, L) belonging to the Laughlin se¬ 
quence [40l 01] , with the same chemical potential /i = 
jiiL = Vhk-F^L = I^R = They are kept at two 

different temperatures Tl > Tr and coupled by a tun¬ 


neling region, as shown in Fig. The two quantum 
Hall systems (QHS) have counterpropagating single edge 
channels with Hamiltonian (in this work h = kB = 1) 

Ha = — [ dxpl{x) = [ dx{dx4>a{x)f . (1) 

Ua J iT^Va J 

Here, Vo, is the propagation velocity of the mode and Va = 
l/rUo,, with rria > 1 an odd integer [40]. The case rria = 1 
corresponds to an integer QHS, while > 1 describes 
fractional quantum Hall liquids. In the second expression 
of Eq. 0 the electron density pa{x) is written in terms 
of the chiral bosonic particle-hole collective excitations 
field (j)a{x). Using bosonization technique [4^ [43] the 
electron operator can be also expressed in terms 

of 

V’a(rc) = ^ (2) 

v27ra 

with a a short distance cut-off and Ta the so-called Klein 
factor [42]. The index a = R(+),I/(—) indicates also the 
direction of propagation and is the associated Fermi 
momentum. We assume that the two QHS are tunnel 
coupled with a tunneling Hamiltonian 

Hj,=aJ dxf{x)i;l,{x)i;L{x) + H.c. , (3) 

where f{x) describes the shape of the tunneling re¬ 
gion [44] and A the amplitude strength [45| [46] . In the 
following we will consider a series of multiple n point-like 
contacts equally-spaced m with distance d [45[ [48ll5Q] , 
with f{x) = S{x—jd)/n. By properly acting on the 

gate voltages of the QPCs one can manipulate and tune 
their transmissions, selectively opening or closing some 
of them. The electric charge Jc and heat Jq currents 
can be written in terms of particle and energy variations 
Jn = {Nr — Nl)/2 and Jr = {Hr — Hr)/2 as 


Jc 


-eJi 


N 


Jq — J\ 


H 


IiJn, 


( 4 ) 


with = j dxpa{x) the particle numbers on each edge. 
The averages (...) are taken over the equilibrium states 
of the left and right QHS with respect to their temper¬ 
atures Tl, Tr EH- They can be computed at lowest 
order in the tunneling, using standard perturbation tech¬ 
niques [5^ l53] 

Jc = —2ie|Ap J dxdx' J drf{x)f{x') sin[2^F(^ — x')] 
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and 
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Here, A = A/( 27 ra) and /cp = (A^F,i? + ^f,l)/ 2 the 
average Fermi momentum. In Eqs. we intro¬ 
duced the function with Wait) = 

(0Q,(t)0Q,(O)) —(0^(0)) the bosonic correlator given bv[42l 

[SI Hi [58] 


Wa{t)=\n 


F 

(1 + ^ 


2 

r2 1 


1(1 + icoct) 


( 7 ) 


with r( 2 ;) the Euler 7 -function and ooc the high energy 
cut-off. It is now useful to introduce the Eourier trans¬ 
form Pm^{E) = f In the energy repre¬ 

sentation we get 


2 

Jc = — J dxdx' J dEf{x)f{x') sin[ 2 fcF(^ ~ x')] 


^2ifkF{x-x') 


PmL{E)PmR{ E) . 

+00 


(8) 


Jq = J dxdx' J dEf{x)f{x') cos[ 2 /cf(x — x')] 

—00 

XP^,{E)P^,{-E) . (9) 




In the scaling limit uJc/Ta ^ I 5 EmdE) can be conve¬ 
niently recast as Pm^{E) = Vm^{E)na{E) with 


^mAE) 


(27r)"^° /F, 

Wcr(TOa) \Wc 


THq, — 1 


^ (2 ^ ^ 2 FT^) 


( 10 ) 

which plays the role of an “effective” tunneling density 
of states (DOS) [56ll5^ and na{E) = -h 1]“^, the 

equilibrium Eermi distribution function at temperature 
T, [60|. Note that in the non-inter act ing/integer case 
{uoi = 1, pyicx = 1) the DOS is constant, as for a nor¬ 
mal Eermi liquid, while for the fractional case {E) is 
energy and temperature dependent taking into account 
the non-Eermi liquid nature of the fractional OHS [56ll59] . 
Using the symmetry properties Vm^{E) = P^^(—P), 
and na{E) + na{—E) = 1 the charge and heat currents 
assume the more compact form 


Jc 

Jq 


TP 

27r 


+00 


E 


g{E + g) 


( 11 ) 


EfriL {E)V 

rriR {E)[nL{E)-nR{E)], 

where we introduced the transmission function 


g{E) = J dxdx' f{x)f{x') cos 


— kY[x — X ) 
ji 


( 12 ) 


In this way Eq. (11) takes an analogous form of the well- 
known Landauer-Biittiker expression for transport im- 
[63] with however a renormalized effective DOS P^^(P). 


The transmission is also sensitive to the shape of the 
tunneling constriction within the g{E) function. Eor a 
periodic array of n QPCs, the case of interest here, the 
modulating function is 

n + 2^7i(n-j)cos(2j7?E///) 
gn{E) = - - - 2 - , (13) 

with the dimensionless quantity g = k^d. 

One immediately recognizes that gi (E) = 1 in the case 
of single QPC, while oscillating functions of the form 
cos{2jgE/g) appear for multiple QPCs n>2. 


III. INTERFERENCE INDUCED 
THERMOELECTRIC PHENOMENA 

A. Thermoelectric switching and heat current 
interference 


We start the discussion with the two QHS at integer 
filling factors ul = ur = 1. Here, Eq. ( 11 ) reduces to 
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(14) 


Explicit calculation, inserting Eq. (13), leads to 
47r|Ape/i 
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where 




TT 


l2{x,y) = y 


cosh(7rx) cosh(7rj/) 


TT 


Mx,y) = 2 


sinm(7rx) sinm(7r^) 
X y 


sinh(7rx) sinh(7r^) 


(17) 

(18) 

(19) 


Before discussing the results, in order to make realistic 
predictions, it is useful to restrict the parameter range 
to a set of experimentally reasonable one. We thus es¬ 
timate ^F ^ with ^ 10 nm a typical magnetic 

length of QHS [64] [65]. The velocity of the edge states 
Va is of the order of ^ lO^m/s. Using these values we 
have /i ^ 10 K. We set the temperature range between 
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20 — 300 mK, typical values in which well-developed frac¬ 
tional quantum Hall plateau were measured [411 [64]. Fur¬ 
thermore, we consider distances between the contacts of 
the order 10 — 300 nm (less than the phase coherence 
length at the considered temperatures [39]), correspond¬ 
ing to a dimensionless parameter range ^ 1 — 30. With 
these parameters the QPCs separation d is never much 
larger than the thermal lengths Lq, = Uq,/Tq, [ 66 | [67] • 

We start now to analyze the charge current. Fig. [^) 
shows Jc as a function of r] for different number of QPCs 
at fixed Tl > Tr. In the case of a single QPC n = 1 
(red/dashed-dotted curve), the charge current is always 
zero, because the energy independence of gi{E) = 1, see 
Eq. ([T^ , does not induce particle-hole symmetry break¬ 
ing [59], so that electrons and holes equally contribute 
to transport leading to no net charge current. On the 
other hand, in the presence of a multiple QPC setup 
(n = 2 , n = 10 in the Figure) quantum interference phe¬ 
nomena are responsible for an energy-dependent trans¬ 
mission function ^^(E), effectively breaking in general 
particle-hole symmetry and leading to a non-vanishing 
charge current. Interestingly enough, the charge current 
exhibits an oscillating behavior, switching between posi¬ 
tive and negative values, with principal zeros at 77 = /c 7 r /2 
{k integer). This result suggests that, despite the ther¬ 
mal gradient is fixed, charge can flow either from hot to 
cold or from cold to hot, depending on the parameters of 
the junction only. To shed light on this result, it is useful 
to rewrite Eq. (14) as 




+ 00 


+ 00 

J dEgn{E + lJ.)nL{E) [1 - nR{E)] 


+ e 


j dEgn{E + g)nL{E) [1 - nR{E)] i , (20) 


with na{E) = na{—E) representing the Fermi distribu¬ 
tion for holes and where the first (second) line represents 
the electron (hole) contribution to the charge transport 
from left to right. In this picture, transport is due ei¬ 
ther to electron tunneling, i.e. nL{E) [1 — nR{E)] ^ 0 , 
or to hole tunneling i.e. ul^E) [1 — nR{E)] ^ 0 . Obvi¬ 
ously, these two cases differ for the sign of the carriers, 
as shown in Eq. (20), so that if particle-hole symmetry 
is present, no charge current is expected. In order for 
the charge current to be finite, the transmission func¬ 
tion must break the particle-hole symmetry, differently 
weighting electron and hole contributions. This cannot 
happen for a single QPC with gi{E) = 1 , as schemati¬ 
cally shown in Fig. § 1 ). However, if n > 1 QPCs are 
present, the energy dependent transmission function can 
promote electron tunneling with respect to hole tunneling 
or viceversa, inducing either negative or positive charge 
current respectively. In particular, transitions between 
positive and negative values of the charge current oc¬ 
cur for T] = /c 7 r/ 2 , where gn{E g) = gn{—E + /i): at 
these specific values particle-hole symmetry is restored. 
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Figure 2 . (Color online) Tunneling charge current Jc in 
units of — 2'K\\\^eg/(jjl for ul — ur — I, as a func¬ 
tion of ?7 = k^d. Panel a): different curves represent differ¬ 
ent numbers of QPCs with temperature Tr = 20 mK and 
Tl = 300 mK. The n — 2 curve displays oscillations of pe¬ 
riod TT, modulated by a non-monotonic envelope function that 
reaches its maximum for g ~ 4 . 57 r. More complicated inter¬ 
ference patterns arise for the n = 10 curve. Note that the 
single QPC (n = I) doesn’t break particle-hole symmetry, re¬ 
sulting in a null charge current. Panel b): Density plot of 
charge current for n = 2 QPCs as a function of g {x axis) 
and Tl/Tr {y axis). Tr is fixed and equal to 20 mK, while 
the ratio Tl/Tr goes from I to 15 (corresponding to a maxi¬ 
mum temperature of 300 mK). The figure shows an increasing 
or decreasing monotonic behavior of the charge current as a 
function of the ratio Tl/Tr at fixed 77 , depending on the sign 
of Jc- The oscillating behavior described in panel a) is also 
visible. In both panels the chemical potential is set to /j = 10 
K. 


i.e. electrons and holes contribute equally giving a null 
charge current signal, as shown in Fig. |^). Quite inter¬ 
estingly the sign of the current switches between positive 
and negative values despite the presence of a fixed ther¬ 
mal gradient direction Tr > Tr. This simply reflects a 
change in the majority of the carriers: either electrons. 
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Figure 3. (Color online) Schematic representation of elec¬ 
tron and hole contributions to the charge current in the 
presence of a thermal gradient (here Tl — 300 mK, Tr = 
20 mK, and ji — 10 K). The electron contribution gn{E + 
/j.)nL{E) [1 — nR{E)] is represented in red, while the hole one 
gn{E + fEjfiL^E) [1 — fiR^E)] is represented in blue (see Eq. 

The insets show the transmission functions gn{E + g). 
(a) Single QPC: the transmission function is energy indepen¬ 
dent and electrons and holes contribute equally, thus giving 
Jc = 0 . (b-d) 2 QPCs. (b) g = Att: despite the transmission 
function is energy dependent, it does not break particle-hole 
symmetry, so that Jc = 0. (c) 77 = dvr — 7 r/ 4 : particle-hole 
symmetry is broken and electron tunneling is enhanced with 
respect to hole tunneling, giving Jc < 0. (d) 77 = dyr + 7 r/ 4 : 
particle-hole symmetry is broken and hole tunneling is en¬ 
hanced with respect to electron tunneling, giving Jc > 0. 



g 

b) 



0.0 27r 47r Ott Stt IOtt 


77 


Fig. §), or holes, Fig. [^). This argument remains valid 
also for 77. > 2, but the presence of higher harmonics shifts 
the position of maximal intensity, as shown in Fig. |^). 

Acting on the parameter 77 one can therefore switch the 
sign (and thus the flow) of the charge current: this can 
be used to implement a device that, exploiting quantum 
interference, allows to selectively switch the charge flow 
induced by a fixed thermal bias. To complete the de¬ 
scription, we present in Fig. [^) the density plot of the 
charge current for 77. = 2 as a function of the tempera¬ 
ture ratio TrITr and 77. The switching behavior of the 
charge current is stable against temperature variations. 
The oscillations as a function of 77 have the same zeros 
also with varying temperature. The interference patterns 
are modulated by an envelope function which moves to¬ 
wards higher 77 values while lowering the temperature ra¬ 
tio. They show a power law behavior {TrITr)^ varying 
the thermal gradient and as a function of 77 a dephas¬ 
ing envelope which scales as I/77 for large 77 values, (see 
Eq. (p!^). Note that the crossing to an exponential de¬ 
phasing dependence would be present only at much larger 
temperatures/QPCs separation, out of the considered pa¬ 
rameter range m- 

The arguments exposed above explain also the oscillat¬ 
ing interference patterns of the heat current Jg shown in 
Fig-i Indeed, as done for the charge current, one can 


Eigure 4 . (Color online) Tunneling heat current Jg in units 
of Jq^ = 27r|A|^/i^/a;c for ul = = 1. The chemical po¬ 

tential is set to /i = 10 K. Panel a): Jg as a function of 
77 = kpd. different curves represent different number of QPCs 
with temperature Tr = 20 mK and Tl = 300 mK. As for the 
charge current, oscillations of period tt are visible in the n = 2 
signal. Here heat current oscillates around a non vanishing 
mean value equals to half the value of the single QPC heat 
current. Note that the modulating function of the oscillating 
part of the signal changes its sign at 77 ~ 4 . Stt, displaying a 
sort of beat. More complicated interference patterns arises for 
the 77 = 10 curve. Panel b): density plot of heat current for 
77 = 2 QPCs as a function of 77 {x axis) and TlITr (y axis). 
Tr is fixed and equal to 20 mK, while the ratio TlITr goes 
from 1 to 15 (corresponding to a maximum temperature of 
300 mK). The heat current is an increasing monotonic func¬ 
tion of the ratio Tl/Tr for all values of 77. The oscillating 
behavior described in panel a) is also visible. 
































































6 


rewrite the heat current as 
(' +00 

27r|Ap 


Jq = 


+ 00 


J dEEgn{E + g)nL{E) [1 - nR{E)] 


J dE{-E)gn{E + fi)nL{E) [1 - nR{E)] \ . (21) 


Contrary to Eq. (20), electron and hole contributions 
add up, because the heat current is insensitive to the 
charge of the carriers. The main features are repre¬ 
sented by the presence of minima and maxima as a func¬ 
tion of 7 ^, around the mean value proportional to Ii in 
Eq. (16). Eor n = 2 they coincide with 77 = {2k + l) 7 r /2 
and T] = kiT respectively, and correspond to values at 
which the transmission function at zero energy has a 
minimum or a maximum respectively. These are pre¬ 
cisely the values that give zero charge current. Then 
the heat current has a maximum or a minimum if, due 
to quantum interference, electrons and holes have high 
or low transmission respectively, as schematically repre¬ 
sented in Eig. Generalizing to n contacts, there are 
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Figure 5. (Color online) Schematic representation of electron 
and hole contributions to the heat current (in units of fi) 
in the presence of a thermal gradient (here Tl = 300 mK, 
Tr — 20 mK, and /x = 10 K). The electron contribution 
Egn{E-\-iJL)nL{E) [1 — nR{E)] is represented in red, while the 
hole one {—E)gn{E -h ii)fiL{E) [1 — fiR{E)] is represented in 
blue (see Eq. ([2T|). Contrary to the charge current, the total 
heat current is given by the sum (and not the difference) of 
electron and hole contributions, see Eq. ( |21[ ). The insets show 
the transmission functions gn{E -h g). (^ Single QPC: the 
transmission function gi{E) — 1, giving the maximum heat 
current, see Fig. [^). (b-d) 2 QPCs. The heat current shows 
a sequence of maxima and minima (compare with Fig. 0 - (b) 
Tj = 27r —7r/2: both electrons and holes have low transmission, 
giving a minimum of Jq. (c) 77 = 27r: both electrons and 
holes have high transmission, giving a maximum of Jq. (d) 
77 = 27r + 7r/2: both electrons and holes have low transmission, 
giving a minimum of Jq. 


multiple paths whose phase differences are always mul¬ 
tiples of 2rj. This explains the more complicated inter¬ 
ference patterns shown in Fig. for 77 . = 10. Similar 


curves are obtained for different values of temperatures 
as shown in the density plot of Fig.j^b) for n = 2. Here, 
the magnitude of the interference patterns increase with 
temperature following the power law {TrITrY increas¬ 
ing the temperature gradient. Note that in both panels 
at 77 bigger than a critical value r]c the interference paths 
have a phase shift of tt, with rjc becoming lower and lower 
by increasing the number of QPCs. 

Mathematically, this is due to a change in the sign 
of the envelope function X 2 in Eq. (19). On a more 
physical ground, one can observe that states contribut¬ 
ing to heat transport are mostly distributed around 
E ^ ±T = ±{Tl + Tr)/ 2, as can be argued from Eig. 
HJa). Therefore, the behavior of the heat current de¬ 
pends on whether these states have high or low transmis¬ 
sion, that is, if gn{E -\- ^ maxima (min¬ 

ima) of the heat current are expected. We focus on 
the simple case 77 = 2 , where the transmission function 
shows an oscillating pattern with period AE = irg/r]. 
Then, consider what happens for values of 77 multiple of 
TT, corresponding to maximum transmission at zero en¬ 
ergy, i.e. g 2 {E ^ = 1. In this case, as long as 

AE ^ 47 rT the transmission function is slowly oscillat¬ 
ing so that states contributing to heat transport have 
high transmission, since g 2 {E + g)\ 


E^T 


1 , and a max¬ 
imum of the heat current appears. On the other hand, if 
AE ^ 47 rT the transmission function rapidly oscillates, 
so that gn{E + can be significantly smaller than 

one, leading to a minimum of the current. This mech¬ 
anism is related to the presence of thermal dephasing 
and induces an exchange between maxima and minima 
by increasing AE, that is, by increasing 77 , as shown in 
Eig. 1^ a). One can roughly estimate the crossover as 
AEc ~ 47rT, which gives 77 c ~ and corresponds in¬ 
deed to a critical length dc ~ Et, with Lr = the 
characteristic thermal length. Note that in the case of 
77 > 2 QPCs the presence of additional modulations in 
the transmission function leads to a decrease of the crit¬ 
ical value 77c. 


We now comment on the tunneling between different 
filling factors er ^ el with fractional QHS. Eor sake of 
simplicity, we focus on the two bars with er = 1 and 
er = 1/3. Generalization to other filling factors [4T1 [54] 
are straightforward. Also in this case one can calculate 
Eq. ( pr] ) in analytical form (not quoted). Eig. [^) shows 
the charge current Jc and Eig.|^) shows the heat current 
Jq as a function of 77 at fixed TrITr = 15 for different 
number of QPCs. The oscillating behaviors are again 
present, refiecting the interference patterns. The charge 
current has the same positions of the zeros of the integer 
case. However, the currents show a faster decreasing of 
the visibility, while increasing 77 . This is an hallmark of 
the fractional nature refiected in the peculiar behavior of 
the effective DOS P^^(E). Indeed, the latter acquires 
an energy and temperature dependence, with the well- 
known power-law dependence of a non-Eermi liquid, thus 
modifying the envelope function. Here also the scaling 
behavior with temperature is modified, which is refiected 
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T] 
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Figure 6 . (Color online) Charge (panel a) and heat (panel b) 
currents for ul = 1 and = 1/3 as a function of 77 = kpd 
at fixed temperatures Tr = 20 mK and Tl = 300 mK. Units 
of Jc or Jq are respectively = 27r|A|^e/i^/ct;cand Jq^ = 
27r|A|^/i^/ct;c. Different curvers refer to different number of 
QPCs. The chemical potential is set to /x = 10 K. 


in a strong suppression of the signal for both Jc and Jq . 


B. Heat current rectification 


We now focus on thermal rectification effects. We con¬ 
sider the heat current Jq in two different (and oppo- 

site) configurations: forward Jq = Jq{Tl = = 

Tcoid) and backward = Jq{Tl = Tcoid,!^ = Thot), 
exchanging the temperatures Thot and Tcoid- One can 
conveniently define a rectification coefficient as pUl l 68 ] 




J. 


if) 


rib) 


( 22 ) 


Recalling Eq. ( pT] ), for constant DOS T^^(T') and 
{E) (integeiyn^on interacting case) one obtains vq = 


1 , where the difference between forward and backward 
heat currents is a simple change of the sign. We now 
show that rectification effects with rg 7 ^ 1 are possible 
in the case of tunneling between a fractional QHS and 
an integer one, achieving also large rectification coeffi¬ 
cients. This phenomenon is related to the energy and 
temperature dependent effective DOS (T) proper of 
a fractional filling. In this case, the difference in the 
DOS between ur and i'l breaks the left-right symmetry, 
see Eq. (11), allowing for rg 7 ^ 1. Eor sake of simplic¬ 
ity, hereafter we restrict the discussion to the tunneling 
between ur = 1/3 and vr = 1. However, provided that 
all results and conclusions remain valid and can 
be easily generalized to other filling factors belonging to 
the Laughlin sequence. We point out that for larger dif¬ 
ferences rriR — rriL one would obtain stronger rectification 
effects {e.g. vr = 1/5-z^l = 1 gives higher rectification 
coefficient compared to the case vr = 1/3-= 1)- 

In the case of single QPC it is possible to get a simple 
expression for the rectification coefficient 


_ 17r^ - lOr^ - 7 
^ 7t4 + 10t2 - 17 ’ 


(23) 


in terms of r = Thot/T/oid- Note that for Thot = T^oid 
one has rg = 1 as expected. On the other hand, in the 
limit Thot ^ T^coid (t ^ 1 ) the rectification coefficient 
saturates to the value rg ^ 17/7 = 2.43. 

Moreover, it is interesting to study possible enhance¬ 
ments of the rectification effects due to quantum inter¬ 
ference phenomena arising in the presence of n multi¬ 
ple QPCs, which, as we have shown in the previous 
Section, have a strong impact on heat transport. In 
Fig#) we show the rectification coefficient as a function 
of Thot/T/oid at fixed Tcoid = 20 mK, for different number 
of QPCs. The case of single contact (solid line in the fig¬ 
ure) is an increasing function of r as reported in Eq. (23). 
The different qualitative behavior of rg in the presence of 
n > 1 QPCs is due to the interplay between interactions 
and interference effects. Remarkably, a pronounced peak 
can be observed, which drifts to lower temperatures, with 
higher values, while increasing the number of contacts n. 
In Eig.l^) we report a density plot of rg as a function of 
T] = k^d {x axis) and Thot/77oid {v axis). Interestingly, 
the optimal condition for the enhancement of vq corre¬ 
sponds to values of 77 for which no charge current flows 
in the system. 

We have shown that efficient heat rectification perfor¬ 
mances can be achieved by increasing the number of 
QPCs. There are however limitations, related to the re¬ 
quirement that phase coherence is preserved throughout 
the tunneling paths, thus giving rise to quantum inter¬ 
ference phenomena. This constrain limits the number of 
QPCs that can be created to n < kn/d^ with the in¬ 
elastic mean free path, that can be of the order of several 
jiim in QHS. 
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IV. CONCLUSIONS 


We studied charge and heat transport in two 
temperature-biased QHS coupled by a tunneling region. 
We showed that, when the tunneling is realized via a 
series of point-like contacts, an interference mechanism 
take place for both charge and heat currents. A mul¬ 
tiple QPCs geometry can effectively break particle-hole 
symmetry, leading to a finite thermoelectric charge cur¬ 
rent. Remarkably its sign, that is, the charge current 
flowing either from hot to cold or viceversa, is governed 
by quantum interference and can be manipulated. Inter¬ 
ference effects affect thermal transport as well, with the 
heat current displaying peculiar oscillations as a func¬ 
tion of the distance between the QPCs. We explained 
both the heat current oscillations and the thermoelectric 
switching in terms of different transmission functions for 
electrons and holes, due to the different tunneling paths. 
Moreover, heat rectification can be achieved when con¬ 
sidering two fractional QHS with different filling factors, 
due to the anomalous non-Fermi liquid tunneling density 
of states of the Laughlin state. Despite heat flow is recti¬ 
fied already in the single contact geometry, the presence 
of multiple QPCs allows to exploit the interferometric 
properties to find an optimal working condition for the 
enhancement of the heat rectification effects. 


Figure 7. (Color online) Panel a): heat rectification coefficient 
rq for multiple QPC geometry as a function of Thot/Fcoid for 
?7 = Stt. Different curves refer to different number of QPCs. In 
the case of single QPC (solid line) rq saturates to 2.43 as one 
can see from Eq. (23). Panel b): density plot of the heat rec¬ 
tification coefficient as a function of rj {x axis) and Thot/Fcoid 
{y axis) for n = 10 QPCs. The ratio Thot/Tcoid goes from 1 
to 15, with fixed Tcoid = 20 mK. Larger rectification effect oc¬ 
curs when ?7 = /cTT and is linked to interference patterns. For 
each value of k the rectification coefficient is a non-monotonic 
function of Thot/Tcoid, displaying a maximum whose position 
changes with k. For y = Stt and Thot/Tcoid ~ 12 one can 
reach a maximum value of rq ~ 9.6. 


V. ACKNOWLEDGMENTS 


We acknowledge the support of the MIUR-FIRB2012 - 
Project HybridNanoDev (Grant No.RBFR1236VV), EU 
FP7/2007-2013 under REA grant agreement no 630925 
- COHEAT, MIUR-FIRB2013 - Project Coca (Grant 
No. RBFR1379UX), and the COST Action MP1209. G. 
D. thanks also CNR-SPIN via Seed Project PGESE003. 


[1] M. Esposito, U. Harbola, and S. Mukamel, Rev. Mod. 
Phys. 81, 1665 (2009). 

[2] M. Campisi, P. Hanggi, and P. Talkner, Rev. Mod. Phys. 
83, 771 (2011). 

[3] E. Giazotto, T. T. Heikkila, A. Luukanen, A. M. Savin, 
and J. P. Pekola, Rev. Mod. Phys. 78, 217 (2006). 

[4] N. Li, J. Ren, L. Wang, G. Zhang, P. Hanggi, and B. Li, 
Rev. Mod. Phys. 84, 1045 (2012). 

[5] D. Segal, Phys. Rev. E 90, 012148 (2014). 


[6] K. Saito and A. Dhar, Phys. Rev. Lett. 99, 180601 
(2007). 

[7] R. Lopez, S.-Y. Hwang, and D. Sanchez, J. Phys.: Gonf. 
Ser. 568, 052016 (2014). 

[8] B. Sothmann, R. Sanchez, and A. N. Jordan, Nanotech¬ 
nology 26, 032001 (2015). 

[9] M. Garrega, P. Solinas, A. Braggio, M. Sassetti, and U. 
Weiss, New. J. Phys. 17, 045030 (2015). 

[10] Y. Dubi and M. Di Ventra, Rev. Mod. Phys. 83, 131 

( 2011 ). 















9 


[11] G. Benenti, G. Gasati, T. Prosen, and K. Saito, 
arXiv: 1311.4430. 

[12] F. Giazotto and M. J. Martmez-Perez, Nature 492, 401 

( 2012 ). 

[13] M. J. Martmez-Perez and F. Giazotto, Nature Gommu- 
nications 5, 3579 (2014). 

[14] M. J. Martmez-Perez, P. Solinas, and F. Giazotto, J. Low 
Temp. Phys. 175, 813 (2014). 

[15] R. W. Simmonds, Nature 492, 358 (2012). 

[16] A. Fornieri, M. J. Martmez-Perez, and F. Giazotto, Appl. 
Phys. Lett. 104, 183108 (2014). 

[17] J. Matthews, F. Battista, D. Sanchez, P. Samuelsson, and 
H. Linke, Phys. Rev. B 90, 165428 (2014). 

[18] T. Ruokola, T. Ojanen, and A.-P. Jauho, Phys. Rev. B 
79, 144306 (2009). 

[19] D. Segal and A. Nitzan, Phys. Rev. Lett. 94, 034301 
(2005). 

[20] S.-Y. Hwang, R. Lopez, M. Lee, and D. Sanchez, Phys. 
Rev. B 90, 115301 (2014). 

[21] D. Sanchez and R. Lopez, Phys. Rev. Lett. 110, 026804 
(2013). 

[22] T. Ruokola and T. Ojanen, Phys. Rev. B 83, 241404 

( 2011 ). 

[23] A. Mari, A. Farace, and V. Giovannetti, arXiv: 1407.8364 
(2014). 

[24] R. Sanchez, B. Sothmann, and A. N. Jordan, Phys. Rev. 
Lett. 114, 146801 (2015). 

[25] R. Bosisio, S. Valentini, F. Mazza, G. Benenti, R. Fazio, 
V. Giovannetti, and F. Taddei, Phys. Rev. B 91, 205420 
(2015). 

[26] F. Mazza, S. Valentini, R. Bosisio, G. Benenti, R. Fazio, 
V. Giovannetti, and F. Taddei, arXiv: 1503.01601 (2015). 

[27] L. Arrachea and E. Fradkin, Phys. Rev. B 84, 235436 

( 2011 ). 

[28] D. G. Rothe, E. M. Hankiewicz, B. Trauzettel, and M. 
Guigou, Phys. Rev. B 86, 165434 (2012). 

[29] R. Takahashi and S. Murakami, Phys. Rev. B 81, 
161302(R) (2010). 

[30] R. Takahashi and S. Murakami, Semicond. Sci. Technol. 
27, 124005 (2012). 

[31] T. Yokoyama and S. Murakami, Physica E 55, 1 (2014). 

[32] Y. Xu, Z. Gan, and S.-G. Zhang, Phys. Rev. Lett. 112 , 
226801 (2014). 

[33] E. Grosfeld and S. Das, Phys. Rev. Lett 102 , 106403 
(2009). 

[34] H. Aita, L. Arrachea, G. Naon, and E. Eradkin, Phys. 
Rev. B 88, 085122 (2013). 

[35] M. Konig, S. Wiedmann, G. Briine, A. Roth, H. Buh- 
mann, L. W. Molenkamp, X.-L. Qi, and S.-G. Zhang, 
Science 318, 766 (2007). 

[36] A. Roth, G. Briine, H. Buhmann, L. W. Molenkamp, J. 
Maciejko, X.-L. Qi, and S.-G. Zhang, Science 325, 294 
(2009). 

[37] E. Olshanetsky, Z. Kvon, G. Gusev, A. Levin, O. 
Raichev, N. Mikhailov, and S. Dvoretsky, Phys. Rev. 
Lett. 114, 126802 (2015). 


[38] A. Kononov, S. Egorov, Z. Kvon, N. Mikhailov, S. 
Dvoretsky, and E. Deviatov, arXiv: 1505.04535. 

[39] G. Granger, J. P. Eisenstein, and J. L. Reno, Phys. Rev. 
Lett. 102 , 086803 (2009). 

[40] R. B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 

[41] D. G. Tsui, Rev. Mod. Phys, 71, 891 (1999). 

[42] T. Giamarchi, Quantum Physics in One Dimension (Ox¬ 
ford University Press, Oxford, 2003). 

[43] E. Miranda, Braz. J. Phys. 33, 3 (2003). 

[44] We also neglect sub leading crossed tunneling process of 

the form with x ^ x . 

[45] G. de G. Ghamon, D. E. Ereed, S. A. Kivelson, S. L. 
Sondhi, and X. G. Wen, Phys. Rev. B 55, 2331 (1997). 

[46] G. Dolcetto, S. Barbarino, D. Eerraro, N. Magnoli, and 
M. Sassetti, Phys. Rev. B 85, 195138 (2012). 

[47] A more general situation, however, does not affect our 
conclusions and physical results. 

[48] E. Dolcini, Phys. Rev. B 83, 165304 (2011). 

[49] D. Eerraro, G. Dolcetto, R. Gitro, F. Romeo, and M. 
Sassetti, Phys. Rev. B 87, 245419 (2013). 

[50] Note that the tunneling properties do not depend on the 
choice of the origin of the coordinate system, but only on 
the relative distance between the QPGs. 

[51] Note that in this work we consider only the electronic 
contribution to thermal transport. We neglect phonon 
contributions since for the range of temperature consid¬ 
ered these are negligible. 

[52] T. Martin, Les Houches Session LXXXI^ ed. Bouchiat et 
al, Elsevier, Amsterdam (2005). 

[53] A. Kamenev, A. Levchenko, Adv. Phys. 58, 197 (2009). 

[54] D. Ferraro, M. Garrega, A. Braggio, and M. Sassetti, New 
J. Phys. 16, 043018 (2014) 

[55] G. Dolcetto, L. Vannucci, A. Braggio, R. Raimondi, and 
M. Sassetti, Phys. Rev. B 90, 165401 (2014). 

[56] M. Sassetti and U. Weiss, Europhys. Lett. 27, 311 (1994). 

[57] A. Furusaki and K.A. Nagaosa, Phys. Rev. Lett. 88, 
226404 (2002). 

[58] J. von Delft and H. Schoeller, Ann. Phys. 7, 225 (1998). 

[59] G. L. Kane and M. P. A. Fisher, Phys. Rev. Lett. 76, 
3192 (1996). 

[60] The Fermi distribution is here defined having the zero 
energy at the chemical potential 

[61] U. Sivan and Y. Imry, Phys. Rev. B 33, 551 (1986). 

[62] M. Biittiker, Y. Imry, R. Landauer, and S. Pinhas, Phys. 
Rev. B 31, 6207 (1985). 

[63] Ya. M. Blanter and M. Biittiker, Physics Reports 336, 
1166 (2000). 

[64] G. Altimiras, H. le Sueur, U. Gennser, A. Anthore, A. 
Gavanna, D. Mailly, and F. Pierre, Phys. Rev. Lett. 109, 
026803 (2012) 

[65] D. Ghevallier, J. Rech, T. Jonckheere, G. Wahl, and T. 
Martin, Phys. Rev. B 82, 155318 (2010). 

[66] K. Le Hur, Phys. Rev. B 65, 233314 (2002). 

[67] P. Virtanen and P. Recher, Phys. Rev. B 83, 115332 

( 2011 ). 

[68] M. J. Martmez-Perez, A. Fornieri, and F. Giazotto, Na¬ 
ture Nanotechnology 10, 303 (2015). 



